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Abstract 

Background: Lead (Pb) exposure during pregnancy may increase the risk of adverse maternal, infant, or childhood 
health outcomes by interfering with hypothalamic-pituitary-adrenal-axis function. We examined relationships 
between maternal blood or bone Pb concentrations and features of diurnal Cortisol profiles in 936 pregnant 
women from Mexico City. 

Methods: From 2007-1 1 we recruited women from hospitals/clinics affiliated with the Mexican Social Security 
System. Pb was measured in blood (BPb) during the second trimester and in mothers' tibia and patella 1 -month 
postpartum. We characterized maternal HPA-axis function using 10 timed salivary Cortisol measurements collected 
over 2-days (mean: 19.7, range: 14-35 weeks gestation). We used linear mixed models to examine the relationship 
between Pb biomarkers and Cortisol area under the curve (AUC), awakening response (CAR), and diurnal slope. 

Results: After adjustment for confounders, women in the highest quintile of BPb concentrations had a reduced 
CAR (Ratio: -13%; Confidence Interval [CI]: -24, 1, p-value for trend < 0.05) compared to women in the lowest 
quintile. Tibia/patella Pb concentrations were not associated with CAR, but diurnal Cortisol slopes were suggestively flatter 
among women in the highest patella Pb quantile compared to women in the lowest quantile (Ratio: 14%; CI: -2, 33). 
BPb and bone Pb concentrations were not associated with Cortisol AUC. 

Conclusions: Concurrent blood Pb levels were associated with Cortisol awakening response in these pregnant women 
and this might explain adverse health outcomes associated with Pb. Further research is needed to confirm these results 
and determine if other environmental chemicals disrupt hypothalamic-pituitary-adrenaLaxis function during pregnancy. 
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Introduction 

The hypothalamic -pituitary- adrenal (HPA) axis is re- 
sponsible for the regulation of Cortisol production in 
humans. Cortisol plays a major role in neurogenesis and 
is critical to brain development, especially in memory for- 
mation, with the hippocampus being the brain region with 
the highest concentration of glucocorticoid receptors 
[1,2]. In humans, Cortisol exhibits a diurnal secretory pat- 
tern, where levels sharply increase in the first hour after 
awakening and slowly decline across the day till evening. 
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HPA-axis function can be measured by examining the di- 
urnal pattern of salivary Cortisol concentrations and 
changes in these patterns are associated with health and 
disease [3-6]. Some studies suggest that alterations in ma- 
ternal HPA-axis function during pregnancy are associated 
with neurobehavioral deficits [7-9]. 

Gestational Pb exposure is a well-known risk factor for 
adverse infant, child, and adult behavioral and cognitive 
development [10,11]. Animal and human studies suggest 
that lead (Pb) exposure may be associated with alterations 
in HPA-axis function, where increased Pb exposure is 
associated with altered adrenocortical response to acute 
stressors [12-16]. Thus, the toxic effects of Pb may be 
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mediated in part through alterations in maternal HPA axis 
function. 

Despite the known role that Pb and Cortisol play in 
brain development, the relationship between Pb and Cor- 
tisol demonstrated in animals, and potential for maternal 
HPA-axis disruptions to impact child development, we 
are not aware of any animal or epidemiological studies 
that have evaluated associations between Pb exposure 
and ambulatory diurnal Cortisol patterns. Thus, we in- 
vestigated the relationship between maternal bone or 
blood Pb levels and salivary Cortisol concentrations dur- 
ing pregnancy in a cohort of Mexico City women to 
determine if cumulative or recent Pb exposure was asso- 
ciated with altered diurnal Cortisol profiles or daily Corti- 
sol production. 

Methods 

Participants 

Participants for this study were enrolled from an ongoing 
prospective birth cohort in Mexico City. Between July 
2007 and February 2011, we invited pregnant women re- 
ceiving health insurance and prenatal care through the 
Mexican Social Security System (IMSS) to participate in 
our study. The IMSS is funded by the federal government, 
employers, and employees to provide health care to 
private-sector employees and low- to middle-class workers 
and their families. 

To be eligible for participation in the study, women 
had to be <20 weeks gestation, >18 years old, free of 
heart or kidney disease, have access to a telephone, plan 
to reside in Mexico City for the next 3 years, not use ste- 
roids (including glucocorticoids) or anti-epilepsy drugs, 
and not consume alcohol on a daily basis. Institutional re- 
view boards at the Harvard School of Public Health, Icahn 
School of Medicine at Mount Sinai, and Mexican National 
Institute of Public Health. All women provided informed 
consent after study protocols were explained to them by 
research staff. 

Blood Pb measurements 

Venous blood was collected in trace metal free tubes 
from women upon study enrollment during the 2nd tri- 
mester, briefly refrigerated at 2-6°C, and subsequently 
frozen at -20°C until analyzed for Pb using graphite fur- 
nace atomic absorption spectrophotometry (Model 3000: 
Perkin Elmer, Wellesley, MA, USA) [17]. Analysis of 
blinded quality control samples provided by the Mater- 
nal and Child Health Bureau and the Wisconsin State 
Laboratory of Hygiene Cooperative Blood Lead Profi- 
ciency Testing Program demonstrated good precision 
and accuracy. The limit of detection for this method is 
<1 ug/dL and the instrument precision is approximately 
1 ug/dL. Blood Pb concentrations were examined as 
continuous variables (ug/dL) and in quintiles. 



Bone Pb measurements 

Tibia (cortical bone) and patella (trabecular bone) bone 
Pb concentration measurements were obtained ~1 month 
postpartum (mean: 34 days, range: 26-55) using a K-shell 
X-ray fluorescence instrument [18]. The patella and mid- 
tibial shaft of each leg were measured for 30 minutes and 
measures from each leg were averaged by the inverse of 
the measurement variance. Patella and tibia bone Pb mea- 
surements are thought to reflect Pb exposures over the 
span of decades, with tibia measurements reflecting longer 
time spans [19]. Bone Pb levels as low as 0.3 ug/g can be 
quantified using K-XRF; however negative values are 
sometimes observed when the true bone Pb level is close 
to 0 [20]. In order to best accommodate these negative 
values, we created five bone Pb concentration categories 
where the reference category was all values <2 ug/g and 
four additional categories were based on quartiles of 
the values > 2 ug/g. Retaining negative values has been 
shown to be the least-biased method to accommodate 
left-censored data, including bone Pb measurements 
[20,21]. The average imprecision (±SD) of patella and tibia 
Pb values were 6.4 (±2.3) and 6.9 (±2.3) ug/g, respectively. 

Salivary Cortisol measurements 

Saliva sampling provides a non-invasive measure of HPA- 
axis function [3]. Between 14 and 35 weeks gestation 
(mean: 19.7, standard deviation [SD]: 2.4 weeks), pregnant 
women were given verbal and written instruction describ- 
ing when and how to provide five saliva samples each day 
over two consecutive days during the week or weekend. 
Before collecting samples, research staff demonstrated 
what times of the day women should collect their samples 
and how they should complete their sample collection 
diary. Women were instructed to provide samples into 
Salicaps (IBL International, Hamburg, Germany) using the 
passive drool technique upon awakening, 45 minutes after 
waking, 4 hours after waking, 10 hours after waking, and 
at bedtime. Participants were instructed not to eat, brush 
their teeth, or drink liquids for at least 15 minutes before 
providing a sample and not drink caffeinated beverages 
before collecting the first two samples. After sample col- 
lection, women recorded the collection time on the tube 
and in their diary, and were asked to refrigerate samples 
until pickup, after which they were stored at -70°C until 
shipment to Dresden, Germany on dry ice. 

Saliva samples were assayed in the same batch in 
duplicate for Cortisol using a chemi-luminescence-assay 
with sensitivity of -0.16 ng/ml (IBL; Hamburg, Germany, 
Clemens Kirschbaum). Control sera covering at least three 
levels of Cortisol were run during each 24-hour time 
period and intra- and interassay coefficients of variation 
were less than 8%. 

Salivary Cortisol concentrations exhibit a diurnal rhythm, 
rising very quickly after waking and then falling throughout 
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the waking hours. We examined total daily Cortisol pro- 
duction and two features of women's diurnal Cortisol pat- 
terns. The area under the curve (AUC) of the observed 
salivary Cortisol levels is an approximation of total daily 
Cortisol production. It was computed using the trapezoidal 
rule and standardized to the mean value that all women in 
the study spent awake each day (~15 hours) [3]. The rapid 
increase in Cortisol concentrations in the morning, known 
as the Cortisol awakening response (CAR) was estimated 
from the change in Cortisol concentrations between the 
first and second saliva samples of each day. The decline in 
concentrations over the course of the day is known as diur- 
nal rhythm and was estimated using the change in salivary 
Cortisol concentrations between the 1st and 5th samples of 
each day. 

Covariates 

We considered covariates that might be associated with 
both maternal Pb exposure and salivary Cortisol concen- 
trations [3,22]. Maternal age, marital status, years of 
education, parity, and smoking status (never, former, and 
current) was collected using standardized questionnaires. 
The week Cortisol sample collection during pregnancy 
was calculated based on date of last menstrual period. 
Weight and height were measured at the baseline study 
visit and used to calculate body mass index (BMI). We 
measured calf-circumference with a tape measure at the 
1 month postpartum visit and adjusted for it in bone Pb 
models. 

We also collected information on stress or depressive 
symptoms, although we did not adjust for these in our 
primary models since this information was collected in 
the 3rd trimester and these symptoms may be a causal 
intermediate between Pb exposure and Cortisol [23]. 
Maternal stress was assessed using the Crisis in Family 
Systems-Revised survey, which inquires about 11 life 
event domains experienced in the last 6 months (finances, 
legal problems, career, home events, relationships, safety 
in the home, safety in the neighborhood, personal medical 
issues, medical issues with others, prejudice, and authority 
figures/institutions) [24]. A summary score (range 0 to 11) 
was calculated by summing the number of domains that 
women endorsed as having at least one negative life event. 
We measured depressive symptoms using the Edinburgh 
Postnatal Depression Scale [25]. 

Statistical analyses 

Given their right-skewed distribution, all Cortisol mea- 
sures were logio-transformed. We used linear mixed 
models with random intercepts for day and participant, 
assuming an unstructured covariance matrix to account 
for the within-woman correlation of Cortisol concentra- 
tions within and across days [26] . To visualize patterns of 
Cortisol response over the course of a day, we modeled 



Cortisol concentrations as a function of hours since waking 
using a five-knot restricted cubic polynomial spline [27]. 
We created smoothed plots according to quantiles of 
blood and bone Pb concentration to visually assess the im- 
pact of Pb on Cortisol rhythms. 

We report the percent difference in Cortisol AUC ac- 
cording to Pb exposure categories. To estimate the 
association between Pb and CAR, we modeled salivary 
Cortisol concentrations at the first and second samples 
as a function of time since waking (hours), Pb levels, 
their interaction, and covariates. We examined the diur- 
nal slope using the first and fifth saliva samples and the 
same approach used for the CAR. We report the percent 
change in Cortisol per hour according to Pb levels and 
the beta coefficient of the ratio of these percent changes 
in Cortisol per hour according to Pb levels. We restricted 
these two analyses to women who reported collecting 
their samples in the a prior defined sampling windows 
since salivary Cortisol concentrations change rapidly in 
the first hour after waking and incorrectly timed samples 
can misclassify Cortisol profiles (see Additional file 1: 
Table SI for a priori defined windows) [28]. 

Secondary analyses 

We conducted a secondary analysis adjusting for 3rd 
trimester maternal stress and depressive symptoms, the 
day of the week (weekday vs. weekend), the time of 
awakening, or number of hours slept on the prior night. 
We also examined diurnal slopes using the 1st, 3rd, 4th, 
and 5th salivary Cortisol concentrations collected each 
day. Finally, we re-conducted our tibia and patella Pb 
analyses after excluding women whose bone Pb mea- 
surements imprecision values were >10 (n = 48-53) and 
15 (n = 1) ug/g, respectively. 

Results 

A total of 1,054 women enrolled in our study during 
their 2nd trimester of pregnancy. Of these, 937 women 
provided at least one day of valid Cortisol samples. We 
excluded one woman who had 4 out of the 10 highest 
Cortisol values in the entire cohort. A total of 9,183 sal- 
iva samples were provided by the remaining 936 women. 
With the exception of former-smoking (33.8 vs. 39.8%) 
and marital status (57.4 vs. 47.5% married), women in- 
cluded in our analyses were similar to the 118 women 
not included. Our primary blood Pb analyses included 
873-918 women with complete data and our bone Pb 
analyses included 568-594 women with complete data. 

On average, women were 28 years of age, modestly ed- 
ucated (mean: 12 years), married (57%), had no children 
or one child (74%), and never smokers (66%) (Table 1). 
A total of 199 women (21.3%) had blood Pb concentra- 
tions greater than 5 ug/dL and 34 (3.6%) had concentra- 
tions greater than 10 ug/dL. Logi 0 -transformed blood Pb 
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Table 1 Univariate description of demographic characteristics, bone/blood Pb concentrations, and Cortisol sample 
characteristics of 936 Mexico City pregnant women (2007-201 1) 



Variable 

Materna Age 

Materna Education 

Materna Weight 

Materna Height 

Materna BMI (kg/m 2 ) 

Weeks Gestation at Cortisol Collection 

Calf Circumference at 1-Month visit (cm) 

Negative Life Events Score 

Edinburgh Postnatal Depression Scale Score 

Maternal 2nd Trimester Blood Pb Concentration (ug/dL) 

Maternal Tibia Pb Concentration (ug/g) 

Maternal Patella Pb Concentration (ug/g) 

Waking Time (SD in minutes) 

Hour of Sleep (hours) 1 

Marital Status 

Married 

Free Union 

Single/Divorced/Separated 
Number of prior children 
1 

2 
3 

4+ 

Smoking Status During Pregnancy 

Never 

Former 

Current 

Weekday Cortisol Collection' 

No 

Yes 



N (%) 



537 (57.4) 
226 (24.1) 
173 (18.5) 

356 (38.0) 
339 (36.2) 
175 (18.7) 
66 (7.1) 

614 (65.6) 
316 (33.8) 
6 (0.6) 

333 (35.6) 
603 (64.4) 



Mean (SD) 

27.8 (5.5) 

1 1 .9 (2.9) 
64.7 (11.0) 

1.6 (0.1) 
26.9 (4.2) 
1 9.7 (2.4) 
33.6 (3.2) 
3.1 (2.1) 
8.4 (5.6) 

3.7 (2.7) 
2.7 (8.4) 
4.6 (8.6) 
0747 (83) 
8.3 (1.6) 

537 (57.4) 
226 (24.1) 
173 (18.5) 

356 (38.0) 
339 (36.2) 
175 (18.7) 
66 (7.1) 

614 (65.6) 
316 (33.8) 
6 (0.6) 

333 (35.6) 
603 (64.4) 



Min, Max 

18.0. 44.0 
0.0, 21.0 

40.0, 108.0 

1 .4, 1 .8 
1 7.4, 44.7 

14.1, 35.1 
25.0, 47.0 
0.0, 1 1 .0 
0.0, 28.0 

0.5, 23 
ND, 20 
ND, 43 
0305-1300 
0.0, 14.5 



ND: Not detectable. 

1-On first day of collection. 



concentrations were modestly correlated with tibia (Pearson 
R = 0.26, p < 0.001) and patella (Pearson R = 0.37, p < 0.001) 
Pb concentrations. Patella and tibia Pb concentrations 
were modestly correlated as well (Pearson R = 0.38, p < 
0.001). 

As expected, salivary Cortisol concentrations exhibited 
a diurnal shape over the course of the day, where they 
were highest in the hour after awakening, then de- 
creased steeply in the next 5 hours after waking with the 
slope flattening across the later 10 hours of the day 
(Additional file 1: Figure SI). 

Salivary Cortisol AUC concentrations did not vary with 
blood or bone Pb concentrations (Table 2), although the 



Cortisol AUC was suggestively higher among women in 
the two highest quantiles of patella Pb. 

Salivary Cortisol concentrations increased 5-7% during 
the first hour of the day among women in the 1st, 2nd, 
and 3rd quintiles of blood Pb concentrations, but de- 
creased 5 and 8% among women in the 4th and 5th 
quintiles, respectively (Table 3 and Figure 1). The CAR 
among women in the fourth (p: -0.04; 95% confidence 
interval [CI]: -0.11, 0.02) and fifth (|3: -0.06; 95% CI: -0.12, 
0.0) quintiles were lower than women in the first 
quintile. There was not a consistent pattern of relation- 
ships between bone Pb concentrations and CAR slopes 
(Figure 2). 
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Table 2 Adjusted change in Cortisol area under the curve (AUC) according to blood, tibia, and patella Pb 
concentrations in Mexico City women (2007-2011) 



Pb biomarker/quantile N 
Blood Pb 1 

Quintile 1 (0-< 1.8ug/dL) 185 

Quintile 2 (1.8- < 24 ug/dL) 183 

Quintile 3 (2.4- < 3.4 ug/dL) 184 

Quintile 4 (3.4- < 5.1 ug/dL) 182 

Quintile 5 (>5.1 ug/dL) 184 

Continuous, per 2.5 ug/dL 918 
Tibia Pb 2 

Quantile 1 (<2 ug/g) 266 

Quantile 2 (2- < 4.3 ug/g) 82 

Quantile 3 (4.3- < 6.7 ug/g) 82 

Quantile 4 (6.7- < 11.1 ug/g) 82 

Quantile 5 (>1 1.1 ug/g) 82 
Patella Pb 2 

Quantile 1 (<2 ug/g) 241 

Quantile 2 (2- < 4.5 ug/g) 86 

Quantile 3 (4.5- < 7.8 ug/g) 87 

Quantile 4 (7.8- < 1 2.7 ug/g) 86 

Quantile 5 (>1 2.7-43.2 ug/g) 87 



GM salivary Cortisol AUC (nmol-hours) % Difference in Cortisol AUC nmol-hours (95% CI) 



82.0 
88.6 
89.8 
88.8 
83.8 



87.8 
83.2 
89.3 
87.5 
93.4 

86.6 
87.6 
81.8 
97.0 
90.3 



Ref 

8 (-1, 18) 

9 (0, 1 9) 
8 (-1, 18) 
2 (-6, 1 2) 
-1 (-4, 1) 

Ref 
-5 (-14, 5) 
2 (-8, 1 3) 

0 (-10, 10) 
6 (-4, 18) 

Ref 

1 (-8, 1 2) 
-6 (-14, 4) 

12 (1, 24) 
4 (-6, 1 6) 



GM: Geometric Mean, AUC: Area Under the Curve. 

1 - Adjusted for maternal age (continuous years), education (continuous years), body mass index (continuous), marital status (married, free union, or single/ 
separated/divorced), gestational age at time of Cortisol collection (continuous weeks), parity (1, 2, 3, or 4+), and maternal smoking status during pregnancy 
(never, former, or current). 

2- Adjusted for above covariates and maternal calf circumference at 1 -month postpartum. 



Diurnal Cortisol slopes were not different across quin- 
tiles of blood Pb concentrations (Table 4 and Figure 1). 
Women in the top quantile of tibia Pb concentrations 
had flatter diurnal slopes (74% decline per 15 hours) 
than women with tibia Pb concentration < 2 ug/g (77% 
decline per 15 hours), but the 95% confidence interval 
(CI) for this difference included the null value (|3: 0.04; 
CI: -0.03, 0.11). A similar pattern was seen when com- 
paring women in the top (-73% decline per 15 hours) 
and bottom (-77% decline per 15 hours) quantiles of pa- 
tella Pb concentrations (|3: 0.06; CI: -0.01, 0.12), but this 
difference did not reach conventional levels of statistical 
significance. 

Secondary analyses 

Adjusting for 3rd trimester stress or depressive symp- 
toms, time of waking, and hours of sleep did not sub- 
stantively change the pattern of most of our results 
(subset of 709 women). The pattern of results for diurnal 
slopes was similar regardless of whether we used 
2 or 5 samples to characterize the diurnal slope. Our re- 
sults were not substantively changed when we excluded 
women with relatively imprecise bone Pb measurements 
(results not shown). 



Discussion 

Among pregnant women from Mexico City, second 
trimester salivary Cortisol concentrations followed the ex- 
pected diurnal pattern, demonstrating that our study 
protocol was successful in capturing Cortisol daily rhythms. 
Characteristics of this diurnal pattern differed according to 
blood and bone Pb concentrations after adjusting for 
demographic and anthropometric characteristics suggest- 
ing that Pb exposure could alter Cortisol rhythms. How- 
ever, associations between bone Pb concentrations and 
diurnal Cortisol rhythms did not reach convention levels of 
statistical significance. 

The use of multiple Pb exposure biomarkers allowed 
us to finely dissect this relationship with respect to dif- 
ferent windows of exposure. The association between Pb 
exposure and salivary Cortisol rhythms depended on the 
Pb exposure measurement matrix. This is intriguing 
given that these two matrices represent different sources 
and timing of Pb exposure. Blood Pb levels are corre- 
lated with the pool of biologically active Pb and repre- 
sent recently ingested or inhaled Pb, as well as Pb 
released from the skeleton [19]. In contrast, bone Pb 
represents past Pb exposures, including those from de- 
cades ago, that have been deposited into hydroxyapatite 
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Table 3 Adjusted change in Cortisol AM rise according to blood, tibia, and patella Pb concentrations in Mexico City 
women (2007-2011 J 1 

Pb biomarker/quantile % Change in Cortisol per hour (95% CI) B Coefficient (95% CI) 
Blood Pb (n = 873) 2 

Quintile 1 (0- < 1 .8 ug/dL) 5 (-5, 1 6) Ref 

Quintile 2 (1.8- < 2.4 ug/dL) 5 (-5, 16) 0.00 (-0.06, 0.06) 

Quintile 3 (2.4- < 3.4 ug/dL) 7 (-4, 18) 0.01 (-0.05, 0.07) 

Quintile 4 (3.4- < 5.1 ug/dL) -5 (-14, 5) -0.04 (-0.1 1, 0.02) 

Quintile 5 (>5.1 ug/dL) -8 (-16, 2) -0.06 (-0.12, 0.00) 

Continuous, per 2.5 ug/dL N/A -0.02 (-0.04, 0.00) 
Tibia Pb (n = 575) 3 

Quantile 1 (<2 ug/g) 4 (-4, 12) Ref 

Quantile 2 (2- < 4.3 ug/g) -8 (-21, 7) -0.05 (-0.12, 0.03) 

Quantile 3 (4.3- < 6.7 ug/g) 7 (-7, 23) 0.02 (-0.05, 0.09) 

Quantile 4 (6.7- < 1 1 .1 ug/g) -9 (-22, 5) -0.05 (-0.1 2, 0.02) 

Quantile 5 (>1 1.1 ug/g) 8 (-7, 25) 0.02 (-0.05, 0.09) 
Patella Pb (n = 568) 3 

Quantile 1 (<2 ug/g) 2 (-6, 10) Ref 

Quantile 2 (2- < 4.5 ug/g) 1 1 (-3, 27) 0.04 (-0.03, 0.1 1 ) 

Quantile 3 (4.5- < 7.8 ug/g) -1 1 (-22, 3) -0.05 (-0.1 2, 0.02) 

Quantile 4 (7.8- < 12.7 ug/g) -2 (-15, 12) -0.01 (-0.08, 0.06) 

Quantile 5 (>1 2.7-43.2 ug/g) 7 (-7, 23) 0.03 (-0.04, 0.10) 

1- Only includes women who collected their first and second Cortisol samples in the correct time window. 

2- Salivary Cortisol concentrations are modeled as a function of time since waking, Pb biomarker, a Pb biomarker x time interaction term, and covariates. 

3- Adjusted for above covariates and maternal calf circumference at 1 -month postpartum. 




Time Since Waking (hours) 

0-<1.9mcg/dL 1 .9 - <2.9 mcg/dl 

2.9 - <3.9 mcg/dL 3.9 - <5.8 mcg/dL 

>=5.8 mcg/dL 



Figure 1 Adjusted geometric mean salivary Cortisol concentrations as a function of time since waking according to 2nd trimester 
blood Pb quintile . 1 -Smoothed function was derived by modeling the time since waking as a restricted cubic polynomial spline with 918 
women's salivary Cortisol concentrations as the outcome. Only samples collected in the correct time windows are included in this plot. The model 
was adjusted for maternal age (continuous years), education (continuous years), body mass index (continuous), marital status (married, free union, 
or single/separated/divorced), gestational age at time of Cortisol collection (continuous weeks), parity (1, 2, 3, or 4+), and maternal smoking status 
during pregnancy (never, former, or current). 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 

Time Since Waking (hours) 

<2.0 mcg/g 2.0 - <4.5 mcg/g 

4.5 - <7.8 mcg/g 7.8 - <12.7 mcg/g 

>=12.7 mcg/g 



Figure 2 Adjusted geometric mean salivary Cortisol concentrations as a function of time since waking according to 1 -month 
postpartum patella Pb quantile . 1 -Smoothed function was derived by modeling the time since waking as a restricted cubic polynomial spline 
with 594 women's salivary Cortisol concentrations as the outcome. Only samples collected in the correct time windows are included in this plot. 
The model was adjusted for maternal age (continuous years), education (continuous years), body mass index (continuous), marital status (married, 
free union, or single/separated/divorced), gestational age at time of Cortisol collection (continuous weeks), parity (1, 2, 3, or 4+), and maternal 
smoking status during pregnancy (never, former, or current). 



Table 4 Adjusted change in Cortisol diurnal slope according to blood, tibia, and patella Pb concentrations in Mexico 
City women from 2007-201 1 1 

Pb biomarker/quantile % Change in Cortisol per 15 hours (95% CI) B Coefficient (95% CI) 

Blood Pb (n = 918) 2 



Quintile 1 (0-< 1.8ug/dL) 


-77 (-79, -75) 


Ref 


Quintile 2 (1.8- < 2.4 ug/dL) 


-77 (-79, -75) 


0.00 (-0.06, 0.05) 


Quintile 3 (2.4- < 3.4 ug/dL) 


-76 (-78, -73) 


0.03 (-0.03, 0.08) 


Quintile 4 (3.4- < 5.1 ug/dL) 


-76 (-78, -73) 


0.02 (-0.03, 0.08) 


Quintile 5 (£5.1 ug/dL) 


-78 (-79, -75) 


-0.01 (-0.06, 0.05) 


Continuous, per 2.5 ug/dL 


N/A 


0.00 (-0.02, 0.02) 


Tibia Pb (n = 594) 3 






Quantile 1 (<2 ug/g) 


-77 (-79, -76) 


Ref 


Quantile 2 (2-<4.3 ug/g) 


-79 (-82, -76) 


-0.05 (-0.11, 0.02) 


Quantile 3 (4.3- < 6.7 ug/g) 


-77 (-79, -73) 


0.00 (-0.07, 0.07) 


Quantile 4 (6.7- < 1 1.1 ug/g) 


-77 (-80, -74) 


-0.01 (-0.08, 0.06) 


Quantile 5 (S1 1.1 ug/g) 


-74 (-78, -71) 


0.04 (-0.03, 0.11) 


Patella Pb (n = 587) 3 






Quantile 1 (<2 ug/g) 


-78 (-79, -76) 


Ref 


Quantile 2 (2-<4.5 ug/g) 


-77 (-80, -74) 


-0.01 (-0.08, 0.05) 


Quantile 3 (4.5- < 7.8 ug/g) 


-77 (-80, -74) 


-0.01 (-0.08, 0.05) 


Quantile 4 (7.8- < 12.7 ug/g) 


-79 (-81, -76) 


-0.04 (-0.11, 0.03) 


Quantile 5 (>1 2.7-43.2 ug/g) 


-73 (-76, -69) 


0.06 (-0.01, 0.12) 



1- Only includes women who collected their first and second Cortisol samples in the correct time window. 

2- Salivary Cortisol concentrations are modeled as a function of time since waking, Pb biomarker, a Pb biomarker x time interaction term, and covariates. 

3- Adjusted for above covariates and maternal calf circumference at 1 -month postpartum. 
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crystal [3]. Consistent with these two matrices represent- 
ing different periods and sources of exposure, we ob- 
served modest correlations between bone and blood Pb 
concentrations. We speculate that past or concurrent Pb 
exposures - represented by bone and blood Pb levels, 
respectively - could have different effects on HPA-axis 
function. 

Specific periods of vulnerability have been identified in 
prior studies of Pb exposure and child cognition [29], 
but little is known about the potential associations 
between Pb exposures at different periods of life and 
HPA-axis function. Pb exposures earlier in pregnancy 
or before pregnancy could impact maternal HPA-axis 
function and future studies should consider additional 
windows of vulnerability (e.g., 1st trimester). Alterna- 
tively, specific pattern of associations could result from 
pregnancy-induced changes in bone Pb mobilization, 
bone turnover, and HPA-axis function [30-32]. 

We are not aware of other studies examining the rela- 
tionship between environmental chemical exposures and 
daily diurnal Cortisol patterns in any human populations. 
Two previous studies have examined the relationship be- 
tween blood/bone Pb concentrations and the response 
to an acute adrenocortical stressor response in children 
and adults [12,13]. Gump et al. found that elevated ges- 
tational and childhood blood Pb concentrations were as- 
sociated with heightened Cortisol responses in children 
following an acute stressor task [12]. Fortin and col- 
leagues found the opposite association in a group of oc- 
cupationally exposed men [13]. In their study, bone and 
blood Pb concentrations were associated with increased 
adrenocorticotropic hormone (ACTH) concentrations 
and decreased Cortisol concentrations before the applica- 
tion of an acute stressor and increased ACTH: Cortisol 
ratios after stressor application, suggesting that Pb expos- 
ure may be associated with adrenal hypo-responsiveness to 
pituitary stimulation. Consistent with this, experimental 
studies in rats suggest that Pb exposure is associated with 
elevated corticosterone concentrations during pregnancy 
[16]. Studies of adult rodents found that Pb exposure was 
associated with elevated ACTH and corticosterone con- 
centrations, as well as increased anxiety-like behaviors [14]. 

The presented results and prior epidemiological stud- 
ies illustrate potentially important differences in HPA- 
axis function response to Pb exposure that may depend 
on the choice of measure for HPA-axis function (e.g., 
acute stressor vs. diurnal rhythm) and timing of Pb ex- 
posure. The latter is important for fetal development 
since disruption of homeostatic Cortisol function during 
the sensitive prenatal period may have profound effects 
later in life since the developing fetus is more sensitive 
to environmental exposures [33]. Changes in maternal 
Cortisol levels during gestation could alter placental ex- 
pression of enzymes involved in glucocorticoid signaling 



and metabolism, ultimately leading to altered neurobe- 
havior [8,9]. 

Both gestational Pb and stress exposure may increase the 
risk of poorer mental and physical development, reduced 
IQ, and delinquent behavior in children [7,10,34, 35]. The 
neurotoxic mechanisms of Pb exposure are diverse and 
may be outcome dependent and the association between 
Pb and HPA-axis function has been relatively unstudied as 
a potential mechanism of Pb toxicity [36]. A prior study 
found that higher maternal 3rd trimester serum Cortisol 
levels were associated with lower child IQ [7]. While this 
finding may seem discordant with our finding of reduced 
CAR among women with elevated blood Pb levels, the 
Cortisol measures in this and the prior study reflect differ- 
ent features of HPA-axis function. The CAR is the change 
in Cortisol over the first hour of the day, while a single 
serum measurement assesses Cortisol concentrations at 
one time during the day. Additional studies will need to 
address the relationships between Pb, diurnal Cortisol 
rhythms, and child neurodevelopment to determine if 
alterations in HPA-axis function are a mediator of Pb- 
induced cognitive deficits. Future studies should also iden- 
tify which patterns of maternal diurnal Cortisol rhythms 
during pregnancy are associated with adverse child 
neurodevelopment. 

While salivary Cortisol concentrations are a non-invasive 
method of measuring HPA-axis function, there are sources 
of variability that may impact the accuracy and precision 
of our results. First, Cortisol exhibits both within-day and 
within-person variability that we attempted to reduce and 
account for the by collecting 10 timed samples across two 
days and using appropriate statistical models. Second, ac- 
curate estimation of the CAR and diurnal slope requires 
that the time of saliva sample collection be reported cor- 
rectly. Seventy-percent of women who reported collecting 
their Cortisol samples in the acceptable time windows had 
at least one day with a positive CAR. This is similar to 
prior studies reporting that 75% of persons have increasing 
Cortisol concentrations after awakening [37]. Incorrectly 
reported times of waking or 1st and 2nd sample collection 
could produce a negative CAR that is not due to biological 
variation, but measurement error. This could systematically 
bias our results if socioeconomic factors associated with Pb 
exposure are also associated with protocol adherence; how- 
ever, we have no reason to suspect this is the case [38]. 

Residual confounding may have produced the pattern 
of results we observed. For instance, women with more 
Pb exposure may live in poorer neighborhoods with 
higher crime, resulting in changes in the pattern of 
salivary Cortisol concentrations [39]. This seems unlikely 
given that the pattern of our results was similar when 
we adjusted for negative life events (including neighbor- 
hood safety) and depressive symptoms. Additional con- 
founding may also arise because pregnancy-induced 
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changes in metabolism may alter Cortisol concentrations 
and increase blood Pb levels by releasing bone Pb stores 
through bone resorption. 

Conclusions 

Cortisol awakening response was associated with in- 
creased blood Pb concentrations among pregnant women 
in this cohort. While many studies of endocrine disruption 
have focused on thyroid and gonadal hormones, the 
association between other environmental chemicals and 
Cortisol homeostasis remains unstudied in humans. The 
importance of glucocorticoids in brain development and 
sensitivity of salivary Cortisol levels to a known toxicant 
illustrates the importance of understanding how environ- 
mental chemicals might impact both Cortisol rhythms and 
subsequent child health and development. 
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